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Recently, the ATLAS and CMS experiments have reported the discovery of a Higgs like resonance 
at the LHC. The next analysis step will include the determination of its spin and CP quantum 
numbers or the form of its interaction Lagrangian channel-by-channel. We show how weak-boson- 
fusion Higgs production and associated ZH production can be used to separate different spin and 
CP states. 
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I. INTRODUCTION 

Recently, ATLAS and CMS have reported the discovery of a Higgs-like [T] resonance with a mass around 
126 GeV [2][3]. Decays to 77, ZZ* , WW* , and recently tt [HE] have been established, with coupling strengths 
consistent with the Standard Model [SHE]- While this 'Higgs discovery' is a great triumph for experimental 
and theoretical high energy physics, the detailed study of this new state will require many years of work at the 
LHC and possibly a future linear collider. 

The ultimate goal of such studies will be a term-by-term and channel-by-channel confirmation of the struc- 
ture of the Higgs Lagrangian from data. If we assume that the observed resonance is indeed responsible for 
electroweak symmetry breaking the leading renormalized operators describing the couplings to all massive 
fermions and gauge bosons are fixed. In an effective theory approach higher-dimensional operators will in- 
duce additional (anomalous) Higgs couplings, which can be constrained from measurements [5]. Some of these 
higher-dimensional operators are responsible for the most striking LHC observables — Higgs production in 
gluon fusion and Higgs decays to photons. Both of these operators have mass dimension six before and five 
after electroweak symmetry breaking, but they are not suppressed by a large mass scale. In the presence of 
Yukawa couplings both, fermion and massive gauge boson loops induce higher-dimensional operators suppressed 
by powers of v }10j . 

This structure of higher-dimensional operators suggests that we should have a careful look at the structure 
of operators responsible for the observation of the new 'Higgs' resonance [HHH]. Only after the operator basis 
is fixed we can determine the corresponding couplings. In that sense the current coupling extractions 6 8 rely 
on strong assumptions about the structure of the Higgs Lagrangian or the Higgs couplings. 

The structure of the Higgs Lagrangian is strongly linked to, but not entirely equivalent to the spin and 
quantum numbers of the heavy resonance. An example is the tensor structure of a coupling of a CP-even 
scalar to two W bosons, which can be constructed using the gauge field or the field strength tensor. Only 
the first can generate a mass for the W boson and unitarize the longitudinal Goldstone boson scattering. 
It is well known how to distinguish between different coupling structures, using the Cabibbo-Maksymowicz- 
Dell'Aquila-Nelson jT3l [14] angles in fully reconstructed X —>■ ZZ decays [5| H5Hl9] , We will show in this paper 
that an equivalent way of determining the Higgs coupling structure can be based on weak-boson-fusion (WBF) 
events [TT1 H"2l |2T)H2"!?] . The two methods nicely complement each other, covering a wide range of production and 
decay channels [23]. In addition, we will show how the angular H — > ZZ analysis can be applied to associated 
ZH production with a fully reconstructed decay H — » bb [24] . 

In contrast to these two complete methods based on a full set of angular correlations, the determination of 
coupling structures for example from X — i 77 events is somewhat limited [53]. The same channel makes it 
unlikely that the set of 'Higgs' measurements is due to a spin-1 resonance, which cannot decay to two neutral 
gauge bosons because of the Landau- Yang theorem Nevertheless, we will see that there is a case for a 

dedicated channel-by-channcl LHC study, carefully keeping track of implicit assumptions for example when 
testing spin-2 models with and without form factor corrections. 

In this paper we develop a complete strategy to extract the structure of the 'Higgs' operators in weak boson 
fusion. We start by constructing the appropriate observables and linking them to the X — > ZZ angles. We 
then define the operator basis best suited for weak boson fusion, which is slightly different from the basis used 
for Higgs decays. In the second part of the paper we study jet-jet and jet- A correlations. We show that the set 
of them is sufficient to probe our set of coupling operators and compare the individual distinguishing power for 
the most popular coupling structures. Finally, we add brief comment on the angular correlation in associated 
ZX production, with a hadronic decay X — > bb. 

II. KINEMATICS 

Before we can develop a detailed strategy on how to extract the form of the Higgs interaction Lagrangian at 
the LHC we need to define our theoretical framework, our set of observables, and their relation for example to 
the well-known H — > ZZ analysis. The latter is usually formulated in terms of a set of angles, where the angle 
between the two Z decay planes is one of the most promising observables. It is proto-typical for linking the 
measurement of quantum numbers or Lagrangian operators to angular correlations, leaving any kind of rate 
measurement for an unbiased determination of the corresponding coupling strength. The same approach we 
suggest for an angular analysis of weak boson fusion Higgs events. As we will see below, weak boson fusion will 
offer one additional handle in form of transverse momentum or energy spectra, but its analysis requires some 
care in its interpretation. 
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Figure 1: Kinematic setup for the angular analysis of H — > ZZ events (left) and Higgs events in WBF production 
(right). All angles are defined in Eq.|2| and Eq.Q. 



A. Flipped Nelson 



The 'traditional' observables to measure the coupling structure of a massive state decaying to two weak 
gauge bosons are the Cabibbo-Maksymowicz-DelPAquila-Nelson [T3JQ3] angles. The kinematics for the decay 
X — » ZZ — > 4£ is illustrated in the left panel of Fig. [T] The four Z decay momenta coming from a heavy 
Higgs-like state X are given by 



Px = Pz e + Pz„ 



Pz, 



Pe- +Pe+ 



PZ U =Pn- +Pu+ 



(1) 



For each of these momenta and the beam direction we define unit three-momenta pi in the X rest frame and in 
the two Z e ^ rest frames. Note that the one of the two Z bosons will be far of its mass shell, i.e. p 2 <C m|, but 
this does not pose a problem for the boost into its reference frame. The set of observable spin and CP angles 
are then defined 



Pe- -PZ U 



Pu~ ■ PZ e 



PZe ■ Pb C 



X 



COS(j) e = (pbearo X PzJ ' (pZ,, X Pe~) 



1 A</> = (p e - X p e+ ) ■ (p^- X 



(2) 



The index at the end of each relation indicates the rest frame in which the angles are defined. In the notation 
in Ref. [TS] this corresponds to 4> e — > <i>i and A0 — > <E>. An important feature is that the reconstruction of the 
angles defined in Eq.Q requires a full reconstruction of the 'Higgs' decay at all stages. It does not require 
both Z bosons to be on-shell as long as we can boost into a well-defined center-of-mass frame of the two decay 
leptons. In spite of the suggestive notation the angles (j> and 9 do not stand for opening angles and not azimuthal 
or polar angles. 



As a first illustration we show the A</> dependences for the process pp — > X —> ZZ —> (e + e~)(/i + /i~) in the 
left panel of Fig. [2] Our hypotheses are the three allowed scalar XZZ couplings structures to mass dimension 
six (or five after symmetry breaking) [TT] [57] and a spin-2 operator [15] • The corresponding operators are 
spelled out in Sec. |III A| and |IIIC| 

For the Standard Model coupling we expect this distribution to have a mild modulation, which would vanish 
for large Higgs masses. In contrast, there are clear modulations in Acf> with a phase shift between CP-even and 
CP-odd dimension-5 operators, which can be easily understood from kinematics [111 I12j . 

The X — > ZZ — > 4£ topology and the weak-boson-fusion 'Higgs' production topology 

5i 92 jih ( x -> d d) (3) 

are linked by a crossing symmetry. The labeling of the incoming and outgoing partons as incoming quarks 
gi,2 and outgoing jets ji t % is only meant to allow for a definition of the angles independently of the partonic 
sub-processes. The 'Higgs' decay products can be d = r, W, Z, 7, depending on the channel we are looking 
at [28H30| . For those observables which require a full momentum reconstruction of pd the list of useful Higgs 
decay channels is reduced. 

Our aim is to generalize the angular basis of Eq.([2| to weak boson fusion, guided by the obvious crossing 
symmetry. When moving one of the final state partons to the initial state we replace time-like Z propagators 
with space-like V = W, Z propagators in the i-channel. In this situation we know that the corresponding Breit 
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Figure 2: Normalized A0 distributions for X — > ZZ events (left), WBF production in the Breit frame (center), and 
WBF production in the laboratory frame (right). We show the SM operator XV^V^ (solid), the CP-even dimension-5 
operator XV^^V^" (dashed), the CP-odd dimension-5 operator XV^ uV^ (dotted), and an example for a spin-2 coupling 
(dashed-dotted). The operators are defined in Sees. IbTA] and [ill C| 



frames are the appropriate reference frames. It is denned as the reference frame where the momentum of the 
t-channel particle V is completely space- like and can be reached through a conventional boost. Similar to the 
X — > ZZ case, a full reconstruction of the X decay is required. With this caveat and fixing the directions of 
all the external momenta as shown in Fig. [I] we define a modified version of the five angles in Eq. ([2| , 



cos 0i = p jl ■ p V2 



COS 2 = P J2 ■ PVi 

Vi Breit 



cos 9* = p v ■ p d 

V 2 Breit 



A" 



cos <j>i = (pv 2 x p d ) ■ (pv 2 x p n ) _ cos A(j> = (p qi x p n ) ■ (p q2 x p h ) . (4) 



Vi Breit 



X 



Again, we define the angle between the two planes of the Breit frames as A(f>. We will see below that this angle 
is similar to the azimuthal difference of two jets in the laboratory frame, A(f>jj, even though in Eq.Q <f) and 9 
do not imply azimuthal or polar angles. It should be noted that using these conventions we can define the angle 
= 2(j>i + A<fi which typically lead to modulations for spin-2 resonances. In the notation of the Ref. [12] this 
corresponds to — > $12. 

Since we cannot measure the charge of the four fermions involved in the WBF topology, we need an alternative 
criterion to map the leptonic observables in Eq.Q to an LHC production process. We break the external fermion 
degeneracy by imposing that each incoming quark q\^. largely keeps its direction when the scattering process 
turns into a tagging jet 2x,2- This defines the weak boson 3-momenta 

PVi.a =Pqi, 2 ~Pn,2 , (5) 

where are the incoming partons moving in the beam direction, p qi 2 = ztebcam- 

Technically, the determination of the incoming particle momenta is an issue at hadron colliders (and 
ill-defined in perturbative QCD). We start by reconstructing the final state and map it onto a set of partonic 
momenta. In the absence of additional jet radiation this defines the center-of-mass system for the incoming 
quarks, related to the laboratory frame by a longitudinal boost. The incoming parton momenta are parameter- 
ized as p qi — (E\,0, 0, Ei) and p q2 — (E%, 0, 0, —E2), which we can invert to express E\ and E2 in terms of the 
summed energy and longitudinal momentum entries of the final state particles. The boost from the laboratory 
frame to the center-of-mass system can be approximated by reconstructing the the events' pseudorapidity from 
the detector geometry and massless calorimeter hits. This should be understood as a prescription to obtain a 
set of well-defined (leading order, parton-like) four-momenta rather than reconstructing the actual initial state. 

Among the angles listed in Eq.Q Acj) does not require a reconstruction of the Higgs candidate. It is 
constructed only from the four external partons, two incoming and two outgoing. Comparing Eq.Q and 
Eq.([2]) we see that it corresponds to the angle between the two Z decay planes in X — > ZZ decays. In the 
central panel of Fig. [2] we show it for illustration, testing the same three spin-0 hypotheses as in the left panel 
and properly defined in Sec. Ill A The Standard Model expectation is again relatively flat, with an additional 



residual dependence introduced the interference of longitudinal and transverse amplitudes and by kinematic 
cuts. The two dimension-5 operators show a distinct modulation, in complete analogy to the corresponding 
measurement in X — > ZZ decays, but with a larger amplitude of the modulation. 
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Figure 3: Normalized cos 9 distributions for X — > ZZ events (left), WBF production in the Breit frame (center), and 
WBF production after analytic continuation (right). Again, we show XV" M V M (solid), XV^V 11 " (dashed), XV ^V^" 



(dotted), and the spin-2 example (dashed-dotted), as defined in Sees. Ill A and III C 



A second key observable in the angular analysis of the X — > ZZ topology is the helicity angle 9 efJ _. In the left 
panel of Fig. [3] we show how it discriminates between a spin-0 and spin-2 resonance |15j . The spin-2 Lagrangian 
is denned in Sec. |III C| For the crossed WBF topology the center panel shows the Breit frame version, defined 
in Eq.Q. It is only of limited use to differentiate between the different hypotheses. 

The translation of the H — > ZZ angles into the Breit frame basis of Eq. Q is not the only possible general- 
ization. To avoid a boost into the Breit frame we can alternatively construct a crossed set of momenta simply 
mapping the pp — > Xjj kinematics to the X — ¥ A£ kinematics. This way we can show the distributions in terms 
of the angles defined in Eq.Q: the WBF topology is related to the decay X — s- ZZ —> 4£ via the crossing 
symmetry p qi — > —p Qi . This crossing we can perform by analytically continuing the WBF initial state momenta 
to final state momenta, rendering py time-like. We illustrate this alternative mapping to space-like momenta 
Pv in the right panel of Fig.[3j including basic acceptance cuts. We see that this observable is more distinctive 
than its Breit frame counterpart, but not as clear as the decay angle correlation. 

The different observations discussed in this section gives us confidence that a detailed study of angular 
correlation in WBF events should serve as a powerful supplement to the study of X — > ZZ decays. 

Additional observables have been suggested for a dedicated study of a spin-2 resonance. One of them is the 
sum of the two Vi^ Breit frame angles [12] . Another is the so-called Gottfried- Jackson angle 9qj [3T| 132], 
It is the angle between the momentum of the resonance X in the laboratory frame and one of the X decay 
products in the X rest frame (somewhat in analogy of the top helicity angle). Compared to the angular basis 
of Eq.Q it is very closely related to 6* , as confirmed in Fig. |4j For spin-0 and spin-1 both distributions show 
a flat behavior, as shown in the Appendix. For spin-2 the two angles are essentially equivalent. It will turn out 
that the angle 9* appears to be more stable with respect to UV completions or the presence of form factors 
appearing in the spin-2 theory, so we will focus on 9* in the following. 



B. Hadron collider observables 

While the kinematic observables introduced in Eq.Q are well suited to capture the properties of a Higgs 
candidate and its coupling to gauge bosons, they are not particularly well suited to cope with QCD effects at 
hadron colliders. The reason is that the reconstruction of the Vi 2 Breit frames assumes a complete reconstruc- 
tion of the hard scattering process. This assumption is not realistic in the presence of QCD jet radiation on 
the one hand and underlying event and pile-up on the other hand. 

As a rough guide line we can look at the opening angle between the two Z decay planes in X — > ZZ decays 
again. For the WBF production process it turns into cos A</> = (p qi x pj t ) ■ (p q2 x Pj 2 ), evaluated in the Higgs 
candidate's rest frame. In the laboratory frame each vector (p q x pj) by definition resides in the azimuthal 
plane. In the Higgs candidate's rest frame this changes, but the impact of a modest transverse boost on the 
large longitudinal parton momenta should be small as well. If we compute in the laboratory frame instead, 
it turns into the azimuthal angle between two vectors, each orthogonal to one of the tagging jet direction. 
This is equivalent to the azimuthal angle between the two tagging jets A.<pjj ■ In the right panel of Fig. [2] we 
show this azimuthal opening angle for the WBF production of four Higgs candidates, with the scalar tree-level 
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Figure 4: Normalized cos#* and cos 9gj distributions for a spin-2 resonance coupling only to weak bosons or coupling 
to quarks as well. We show the X decay to tau leptons and to photons. 



and dimension-5 XVV couplings or a spin-2 coupling. As expected, the Standard Model operator predicts 
an essentially flat behavior, while the other two show a distinctly different modulation [TTJ [2T], arguably the 
strongest of the three approaches compared in this figure. 

Encouraged by this observation we define a new set of (essentially) angular correlations for the WBF Higgs 
production process 

qi q-2 -> hh {X -> dd) . (6) 

It simply uses the rapidity differences and azimuthal opening angles between the tagging jets, the Higgs can- 
didate (if reconstructed), and its decay products: 

{Ar) mn , A4> mn } for m,n = ji j2 ,X,d, d. (7) 

The two tagging jets j\ )2 now have to be ordered, so we defined j\ to be harder of the two. Should there be 
additional jets in the event we choose the hardest two jets which qualify as tagging jets following the basic cuts 
defined in Sec. HVl 

Obviously, this set of observables over-constrains the kinematics in the laboratory frame. The dimensionality 
of the phase space is eight, which includes only five independent angles. While the collider observables of Eq.([7]) 
might not be the most obvious choice for testing spin and CP properties of a heavy resonance X, the case of 
A.(f>jj shows their potential, and they are well-defined experimental quantities at hadron colliders. 

One additional observable of interest is the boost which relates the lab frame to the 'Higgs' and the partonic 
Xjj center-of-mass frames. It is related to the kind of parton in the initial state, distinguishing valence quarks 
from sea quarks or gluons. While it is possible to for example construct asymmetries in the two frames, we do 
not find clear links to the coupling structure of a massive particle in weak boson fusion. 



Finally, at least for a cross check it might be useful to add information which is not encoded in the rapidity 
differences or angles in Eq.|7|. One such observable is the transverse momentum of the tagging jet, where 
the splitting probability Pt.l{%) is known to look different for transverse and longitudinally polarized gauge 
bosons [H3 |33] 



P T (x,p T ) oc 



l + (l-x) 2 V \ 1 + (1 



2x {p\ + (1 — a;) m^) 2 2x p\ 



1 



Looking at large transverse momenta px 3> my the radiation of longitudinal V bosons falls off sharper than the 
radiation of transverse V bosons. Even more distinctively, for spin-1 or spin-2 particle production we know that 
these states do not unitarize the partonic processes. This difference can have a huge effect on the transverse 
momenta of the tagging jets [3H[3S], depending on the additional unitarization procedure. 
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Figure 5: Contamination of WBF Higgs events with gluon-fusion events as a function of m™ m . Note that the unphys- 
ical separation of the partonic processes contributing to the order a 2 a g% gg is only possible at leading order and scale 
dependent. The NLO rates are based on Vbfnlo [201136] . 



C. Contaminating sub-processes 

The main advantage of the leptonic decay H — > ZZ — > A£ is that after some basic kinematic cuts the 
backgrounds are negligible. This holds for non-Higgs events, for example from continuum ZZ production, but 
also for signal events which might have different topology or different intermediate states.* 

Once we switch to WBF Higgs production the non-Higgs backgrounds are still small, S/B — 0(1), but not 
entirely negligible. The angular correlations for example of the Z — > tt background process is well understood, 
which means that they can be taken into account in the final analysis Additional partonic sub-process 

contributing to WBF Higgs production are 

qq -> qqH, ggH qg -> qgH qg -> qgH gg -> qqH, ggH , (9) 

all mediated by the effective H gg coupling, i.e. with cross sections of the order a oc o? s g\ gg . Per se, these rates 
are not small. After the basic acceptance cuts, shown in Sec. |IV[ they can exceed the WBF signal. However, 
to remove the QCD-induced Zjj background at order aa 2 s we need to apply a stiff cut on the invariant mass 
of the two tagging jets rrijj > 600 — 800 GeV. This cut reduces the contamination of gluon-fusion Higgs events 
in the WBF sample to a small perturbation. 

In Fig. [5] we show the relative size of the Higgs production rates in gluon fusion and in weak boson fusion. 
Indeed, for m,™ m ~ 250 the two production modes still contribute with similar rates. When we increasing 
to™ 111 towards 1 TeV the combined gluon fusion contribution rapidly falls off. Splitting it into the different 
contributions of order a 2 s gjj gg , as given by the partonic sub-processes, we see that the originally dominant gg 
initial state drops below 20% of the WBF rate around w™ m ~ 500 GeV. The mixed quark-gluon initial state 
is the most dangerous, because it is almost as large as the gluon fusion rate, but it drops only very slowly with 
a remaining 30% correction to the weak boson fusion rate for realistic wi™ m values. In contrast, the purely 
quark-induced sub-process gives a constant correction around 10% and cannot be reduced choosing larger w™ ln 
values. 

In combination, after applying tagging jet acceptance cuts and the usual ti™ 111 cut we roughly expect one 
third of the Higgs signal events to arise from gluon fusion. However, in a realistic analysis we need a minijet 
veto to further reduce the it and QCD Z+jets backgrounds. The original (still largely valid) estimates for the 
survival probabilities against such a veto range around 90% for the WBF signal and 30% for the QCD Z+jets 
background [331 [53] . Gluon fusion Higgs production will be even more strongly suppressed, with a recently 
estimated jet veto survival probability below 20% [371 EH] - This means that after applying such a jet veto the 
gluon fusion contamination will drop from 50% to a negligible 10%. For rrijj > 1 TeV it ranges around 3%. The 



* A potentially dangerous contribution would be from Higgs decays to Zy, where for a dimension-5 HZZ coupling we would not 
expect any hierarchy between the effective gnzz an d gHZ-y couplings. 
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Figure 6: Normalized A(f>jj distributions for Higgs production in gluon fusion, separated by partonic channels to leading 
order. The CP even Higgs is represented on the top figures and the CP odd on the bottom ones. The m™ m cut is 
increased from 200 GeV (left) to 600 GeV (central) and 1 TeV (right). 



interference between gluon fusion and WBF diagrams can be safely neglected, essentially independently of the 
cuts [39 . In our detailed analysis we therefore omit all gluon-induced operators and show the corresponding 
distributions only in the Appendix. 

In addition to the fact that the WBF Higgs event sample will only have a small contribution from gluon 
fusion, there is no problem in carrying a percentage of gluon fusion event through the Higgs coupling analysis. 
For example, for a scalar Higgs boson the CP quantum number entirely fixes the form of the effective Hgg 
coupling structure. As an example, in Fig. [6] we show the A0jj distribution for gluon fusion Higgs production, 
split into the three partonic sub-processes listed in Eq. ([£]). The partonic sub-processes show a very similar 
behavior, essentially independent of the detailed cuts. Once we require at least rnjj n — 600 GeV the curves are 
completely degenerate. The effective Higgs-gluon coupling in the Standard Model has the same structure as the 
CP-even dimension-5 operator for weak gauge bosons, HV flv V fl1 ' , which explains the quantitative agreement of 
these results with the right panel of Figl2l The rn™ ln cut forces all partonic sub-channels into similar kinematic 
configurations, while the Higgs production is governed by the same coupling to mostly t-channel gluons. If there 
should be a good reason to include more gluon-fusion Higgs events in this coupling analysis we can expect stable 
results in the combination of gluon fusion and weak boson fusion, without any complications from different 
partonic sub-processes governed by the parton densities. 

Finally, an additional complication in WBF is the fact that we cannot determine the charge of the incoming 
and outgoing fermions, i.e. we cannot distinguish neutral from charged current contributions. As long as we 
stick to S , [/(2)-invariant operators this does not pose a problem, but there exist Higgs-candidate oper ators 



which are forbidden by the Bose symmetry between two Z bosons and allowed for a W + W~ pair. In Sec. IIIB 
we will come across them. While we can distinguish photon-induced contributions based on the transverse 
momentum spectrum of the tagging jets, WW and ZZ induced WBF events are indistinguishable, with a 
typical relative size 1 : 4, dominated by WW fusion. 



III. LAGRANGIAN 



Measuring the Higgs coupling structure in weak boson fusion is not a theory-independent measurement. 
Instead, we have to implement all different hypotheses in a Monte Carlo event generator, i.e. Madgraph (40j . 
The different predictions then have to be compared to ATLAS and CMS data. In Section |IV| we will limit 
ourselves to signal-only results, noting that in WBF Higgs searches we expect a favorable signal-to-background 
ratio. Moreover, we will only show a limited number of particularly promising distributions to show that we can 
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indeed distinguish all available scenarios for a heavy resonance X. A complete list of all simulated distributions 
can be found in the Appendix. 

A. Spin zero 

The different hypotheses for the nature of a heavy resonance X can be easily ordered by its spin. First, 
we parameterize the relevant interactions of a CP-even or CP-odd scalar X = H,A with gauge bosons in an 
effective, gauge- invariant Lagrangian [TT1 12"71 HT] . 

JO) JO) (0) (0) 

C = g { ° } HV^V - ^ H V^W" - ^ A V^V^ - ^ H G^G^ -\ A G^ v , (10) 

where V — W,Z are the electroweak gauge bosons and G is the gluon. The coupling to the photon can be 
described in complete analogy, but it will not be relevant for our analysis of LHC production channels. For all 
gauge fields the corresponding V*" indicates the field strength tensor and V^ u is its dual. Some of the couplings 
g^ defined this way have mass dimension, i.e. we do not distinguish between order-one coupling factors and 
mass suppression. This simplified notation is appropriate, since throughout our analysis we assume equal rates 
for all coupling operators. In the Standard Model the couplings are = 2My/v and gf" 1 = a s gHtt/3irv 
(in the heavy top limit) [10] . The additional higher-dimensional coupling g% is finite but negligible in the 
Standard Model. 

In an effective theory approach the couplings g^ — 35° \ which carry mass dimension 1/A, can be generated 
through loops. We omit the additional factors 1/A throughout our discussion because it complicates the 
notation and will be of no relevance to our discussion of angular correlations. Nevertheless, we expect a 
hierarchy in the size of scalar 'Higgs' couplings. In general, the far off-shell processes H — > ZZ — > 4£ and 
H — > WW — > l + l~vv come with an additional suppression factor from the leptonic branching ratios. In 
contrast, there is no such suppression for the loop-induced decay H — > 77. In absence of the Standard Model 
coupling g^ the LHC measurements require g^' 22 ^ ~ g^' WW ^ ss lOOOg^ ' 77 ^- The non-diagonal coupling 
^(0,27) j g similarly suppressed in the absence of an on-shell H — > Z7 signal. However, as long as anomaly-free 
SU(2) multiplets run inside the loops all induced couplings tend to be of similar size. This implies that a 

substantial fine tuning would be needed to accommodate a g% -based explanation of the ATLAS and CMS 
results. 

While such an argument based on relative rates is convincing in an effective theory framework with an 
essentially known particle content, the analysis of angular correlations is model independent. It does not 
assume any perturbative field theory and hence applies to all possible fundamental explanations of the 'Higgs' 
sector of the Standard Model. The list of scalar hypotheses which we test later is summarized in Tab. [TJ 

B. Spin one 

The general Lagrangian describing the interaction between a neutral CP-even or CP-odd gauge boson X = 
y( e ) ; anc i the electroweak gauge bosons is 

A = ig[ 1) (W+ W~" - W- v W + ») Y^> + ig^W+W- Y^ 

4-> 



+ g { s ) e^(W+ d P W~)Y^ +tg i 4 L) W+ fl W-^Y^ 



^W+W-^Y^ + efY^) + ig^W+W-Y^" + ig^W+W'^Y^ 



+ g^^°Y^Z v {d p Z a ) + gWytofrZ^Z" . (11) 

It includes g± — g^ describing the coupling to W bosons [S7J and the reduced set </g — for the Z 
couplings : 42J , constrained by the Landau- Yang theorem [55] . All terms in this Lagrangian have mass dimension 

four, except for g± and g^ . Compared to the usual X — » ZZ analysis this WBF Lagrangian is much more 
complex, so we limit our number of hypotheses to one coupling each for CP-even and CP-odd states to W and 
Z bosons, 

l£: si'Vo l~ z : ff^Vo (12) 
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Table I: Properties of the heavy resonance models discussed in this paper. The structure of the coupling to gluons we 
only list for completeness reasons, in the analysis we do not have to take these operators into account. Instead, we show 
them in the Appendix. 



The interactions of with W bosons has the same form as for a Kaluza-Klcin Z boson, i.e. it is a re-scaled 
version of the WWZ coupling in the Standard Model. For the CP-odd F (o) we resort to a special combination 
of dimension-4 operators only. The interactions with Z bosons is the usual complete basis [T4HT6] . These four 
scenarios are summarized in Tab. [ij We will see that the angular correlations in WBF production of a 1^ and 
1^ can be different. The reason is that in our choice of operators the two are at best partially identical, because 
the list of allowed coupling structures of a neutral resonance to the neutral Z is considerably more constrained. 

When extending our analysis of the WBF production process to include decays of the heavy neutral resonance 
for example into a pair of r leptons we need to assume a structure of its coupling to fermions. For simplicity, 
we assume the same vector and axial- vector structure as the Z boson, 

9v = ~ +2sitl 2 8 w g A = ~, (13) 

again re-scaled by a constant factor. Moreover, we only couple the new vector to the third generation, to avoid 
modifications of the WBF production process. 



C. Spin two 

There is a number of commonly used ways to parameterize the interaction of a spin-2 boson |15) . As usually, 
we choose graviton-inspired scenarios to contrast with the scalar Higgs boson. More precisely, we consider a 
massive graviton resonance, which couples to the SM particles through their energy-momentum tensors, 



C 2 



(2) r rpy-V _ (2) r rp, 

3l ^iiv±y 32 LT /ii' J ( 



(2) n rpllV 

33 ^ilv± f 



(14) 



where G pv is the spin-2 resonance and Ty V G . is the energy-momentum tensor of the electroweak gauge bosons, 
the gluon, and the different fermions. For instance, the tensor for the Z boson is given by [13] 



1 Z 



fz P zp 



m%Z»Z v 



with = d^Z,, —d v Z li + ig w cos 9 W 



w+w-) 



(15) 



While conventional graviton excitations have universal couplings strengths 1/A to fermions and gauge fields, 
where A is the scale parameter of the theory, for our test we define four scenarios motivated on the WBF 
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(16) 

A universally-coupled spin-2 particle is heavily constrained by Tevatron data |44j . but we ignore these con- 
straints in our discussion limited to the Higgs sector. All hypotheses shown in Tab. [I] are not meant as a valid 
ultraviolet completion of the Standard Model without a Higgs boson. All they need to serve as are general 
alternative scenarios for the interpretation of a Higgs-like resonance in weak boson fusion. Similar to the spin-1 
case, the couplings to leptons are assumed to include decays of the spin-2 resonance even for 2^ w and 2q CD . 



IV. ANALYSIS 



The study presented in this section is not meant to produce an optimized Higgs couplings analysis using the 
weak-boson-fusion production mechanism. Instead, we study the behavior of the kinematic variables described 
in Sec. |IIB for the different 'Higgs' candidates. For this purpose, we identify the most promising variables 



to separate the different Higgs spins and couplings defined in Sees. Ill AflHI C| Based on this most promising 
subset of observables we will show how the different hypotheses on the nature and the couplings of the heavy 
state X can be distinguished. In Sec. |IVE] we apply a more sophisticated statistical analysis to this problem, 
but eventually we will leave it to the LHC experiments to determine the quantitative impact of our approach. 



For our numerical analysis we implement the interactions Lagrangians Eqs.(10) 
Rules [ 



(111 and fl4| ) into Feyn- 

which provides a UFO model file that we implement into MadGraph5 [4"0], All simulations 



for this section are done at parton level for an LHC energy of 14 TeV. To study the impact of the acceptance 
cuts on our observables we choose three possible set of acceptance cuts: 



1. We start with a minimal set of cuts applied on the outgoing partons 



PT,j 



> 20 GeV 



ARjj > 0.6 



M<5. 



(17) 



2. To suppress backgrounds and a possible contamination with gluon-fusion events we also include 



mjj > 600 GeV 



(18) 



3. As is common in WBF analyses we require a rapidity gap between the two tagging jets 

Arjjj > 4.2 r) h x r] h < . 



(19) 



Note that while we simulate the decay X — > tt [TTJ [35] as a toy channel for including the X momentum and 
the momentum of the decay product we do not apply any cuts on the decay products. In Tab. [H] we present 
the signal and background ratio suppression after cutting on rrijj, Eq.(18l, and after requiring the rapidity 
gap between the two tagging jets, Eq.(19). The signal and the backgrounds are normalized to the rates after 
the acceptance cuts of Eq.(17|. We see that the background suppression and the reduced contamination from 
gluon-fusion are delivered by the rrijj cut. In the following we will advocate that the cuts from Eq.(19) should 



cuts 


signal 

n+ n+ n~ n+ n~ i+ i _ i+ 1 _ 9+ 

U SM U D5 U D5 u D5g u D5g L W L W L Z 1 Z L 


background rrjj 
EW QCD 


rrijj in Eq.(|18[) 
Ar]jj in Eq.(|19|) 


0.442 0.241 0.290 0.072 0.070 0.318 0.420 0.291 0.419 0.161 
0.389 0.068 0.108 0.053 0.053 0.189 0.242 0.134 0.220 0.096 


0.211 0.027 
0.124 0.014 



Table II: Signal and background cross sections for the different 'Higgs' coupling operators, normalized to the acceptance 
cuts Eq.fT7). 
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be avoided because they reduce the number of observables available for the analysis of the Higgs coupling 
structure. For the remaining part of this paper we will only apply Eqs.(17l and (18). 



This reduced set of kinematic cuts has an effect on the distinguishing power of the observable A<j>jj ■ Forcing 
the two tagging jet into a more collinear configuration enhances the modulation for example distinguishing 
the three scalar coupling structures [TTJ [20] • However, we will see that Arjjj itself is a powerful observable 



distinguishing coupling structures, 
Arjjj analysis. 



so this effect can better be exploited through a two-dimensional A(j)jj vs 



A. Tagging jet kinematics 

From the WBF analyses [551 we know what the kinematic features of the two tagging jets and the 
heavy SM Higgs boson are: both tagging jets are distinctively forward and there is only little color activity 
between them [37] , except for final state radiation off the centrally produced Higgs boson. All of this we show 
in Fig. [7 The corresponding distributions for the coupling structures 0q 5 and 0p 5 are very similar, with a 
slight difference due to the absence of the longitudinal amplitude for the 0~ coupling. In the spin-2 case with 
couplings only to weak gauge bosons, 2^ w , the tagging jets are more central. Also including graviton emission 
from quark lines, 2g W+q , shifts the tagging jets into the forward region. This effect is a first indication that for 
a spin-2 interpretation of the Higgs-like resonance the details of the underlying model have a significant effect 
on essentially all observables. A complete set of reference distributions is shown in Fig. [14] in the Appendix. 

Another serious issue arises when we compare the transverse momentum distributions of the tagging jets 
(or the recoiling heavy resonance). The momentum dependence of the spin-2 coupling generates pr.j spectra 
extending beyond a TeV. However, consistent models of a spin-2 resonance will include an additional form factor 
to cut off this tail and render all px distributions more similar to the scalar case [35j . Including distributions 
with an energy dimension, like the transverse momentum of the tagging jets and the central resonance will only 
test different implementations of such a form factor which cannot be derived from first principles. 

The same behavior can be observed in WBF production of heavy W bosons, for example in little Higgs 
models [31]. Just introducing additional heavy gauge bosons appears to predict very large momenta of the 
tagging jets. However, for heavy gauge bosons we can also check a consistent theory including heavy fermion 
partners and indeed find prj distributions very similar to the known Higgs case. To test the stability of or 
results with respect to the UV completion of the spin-2 model we also show results requiring the tagging jets 
to fulfill 

prj < PT X = 100 GcV . (20) 



The starting point of our analysis of 'Higgs' coupling structures in weak boson fusion is the azimuthal angle 
between the two tagging jets [TT]. Its relation to alternative observables we discuss in Sec. II B and specifically 
in Fig. [2] indicating that it should be well suited to distinguish different spin-0 and spin-2 hypotheses. 
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Figure 7: Rapidity of the resonance X (left), rapidity of the leading tagging jet (center), and transverse momentum of 
the leading tagging jet (right) for the different spin-0 and spin-2 hypotheses after the cuts of Eqs.(|17[) and (18 1. The 



hypotheses 0g M (red), 
Tab. [I] 



0p 5 (blue dashed), D5 (green dashed), 2^ w (black), and 2 



EW+q 



(black dashed) are defined in 
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Figure 8: Normalized correlations between the two tagging jets in WBF production of a heavy resonance X. We show 
the difference in the azimuthal angle A(j>jj and the rapidity difference Ar/jj. The spin-0, spin-1, spin-2 interaction 
structures are described in Tab. Nl 0g M (red), 0j 5 (blue dashed), 0p 5 (green dotted); \\ (blue), 1^, (blue dashed), 1^ 
(green), (green dashed); 2g W (black), 2^ w+q (black dashed). For the spin-2 case we also show the distributions after 
requiring p™ ax = 100 GeV (right panels). 



In the upper panels of Fig. M we show the Atpjj distributions for the nine different couplings listed in Tab.|T] 
There are essentially four different patterns in the left three panels: a flat Acj)jj behavior (0g M , l^,2g W ), a 
back-to-back peak at Acfijj ~ tt (0p 5 ,2g W+q ), a preferred angle around tt/2 (0q 5 , 1^), and preferably aligned 
tagging jets A<pjj ~ (1^ z ). Due to the absence of the longitudinal amplitude for the 0p 5 operator we know 
that its distribution follows a 1 — cos 2A(f>jj modulation. On the other hand, 0p 5 follows a cos Atfijj shape from 
the interference between the transverse and longitudinal amplitudes |12j . 

The issue with the proper definition of the spin-2 couplings also appears in this azimuthal observable. First, 
the distributions without (2^ w ) and with (2g W+q ) a coupling to quarks have little in common. Secondly, 



removing the high-energy tale of the jet momenta requiring Eq.(20) strongly affects the full spin-2 distribution 
2g W+q , changing it from distinctively peaked at A(j>jj — tt to essentially identical to the Standard Model Higgs 

scalar 0g M . This is a general feature which we can trace back to the helicities probed by the spin-2 amplitude: 
once we limit it to relatively small energies of external particles it probes exactly the same helicity structure 
as the spin-0 Standard Model operator [12"] . 

In the lower panels of Fig. [8] we show the rapidity difference Ar/jj between the two tagging jets. In the 
standard WBF analyses we require Ar/jj > 4.2, cutting away most of the spin-1 and spin-2 events and keeping 
mostly the Standard Model Higgs 0g M events. Skipping this cut we see that there are three distinct groups 
of curves, with maxima around Ar/jj ~ 2 (2g W ), Ar/jj ~ 4 (0 D5 , l v , 2g W+q ), and Ar/jj ~ 5.5 (0g M ). Asking 



for tagging jets with limited energy, Eq.( |20[ ), turns both spin-2 distributions into an exact copy of the 0g M 
predictions. In this case only two distinctly different patterns survive, SM-Higgs-like and slightly less forward. 

These rapidity differences are strongly correlated with the tagging jet rapidities shows in Fig. [7j Only very 
little information remains in the additional observation of rjj^i once we exploit Ar/jj. This is an effect of the 
WBF ansatz we are using for all resonances hypotheses — if any of the processes under consideration included 
a mixed quark-gluon initial state the boost from the laboratory frame to the center-of-mass frame of the two 
tagging jets might well turn out useful. 

The main advantage of the two observables discussed in this section is that they do not require the recon- 
struction of the heavy resonance. In Sec. IV E we will compare the two observables A<f>jj and Ar/jj in their 



distinguishing power testing some spin-0 and spin-2 coupling structures. A complete set of distributions can be 
found in Fig.[l5]in the Appendix. As expected, the study of decay independent jet correlations is a very useful 
starting point for a 'Higgs' coupling analysis |llj . However, the predictions from a spin-2 model have to be 
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taken with a grain of salt, because they are not stable with respect to (slight) modifications of the underlying 
model. 



B. Higgs-jet correlations 



Adding information on the 'Higgs' momentum requires the reconstruction of the heavy resonance. This means 
that the most promising channel for such an analysis would be the decay X —¥ 77, possible supplemented by the 
approximate reconstruction of X — > tt decays. The latter requires additional assumptions about the couplings 
of the heavy resonance to fermions. This is most obvious in the spin-1 case, where in Eq.(13) we assume that 



the heavy neutral resonance couples to the third generation in complete analogy to a Standard Model Z boson. 

This kind of assumption has for a long time been a problem in the interpretation of Z' searches, so it would 
be beneficial if we could apply an analysis of the coupling structure without considering fermionic couplings. 
Similarly, for WBF Higgs production the spin-2 states coupling to the energy-momentum tensor makes specific 
assumptions about the production vertex and about the decay vertex. It would be preferable if these two 
assumptions could be kept separate. 

A way of including significantly more information than just the jet-jet correlations discussed in the previous 
section and at the same time making minimal assumptions about the decay of the heavy resonance is to study 
angular correlations between the tagging jets and the reconstructed resonance X. In this second step of the 
angular analysis we will focus on the observables Arjjx and A(f)jx, where j denotes the tagging jet with the 
larger transverse momentum. This ordering of the tagging jets induces a difference compared to the Breit-frame 
analysis described in Sec II A[ where the tagging jets are ordered by rapidity. 



In the upper panels of Figj9] we show the rapidity difference between the reconstructed heavy resonance and 
the leading tagging jet. As expected from the tagging jet distributions shown in Sec. |IV A| thc Standard Model 
Higgs 0g M is peaked at large differences Arjjx ~ 2. The unitarized spin-2 resonance has a characteristically 
flat behavior. The Arjjx distributions for all spin-1 coupling structures as well as the dimension-5 scalars are 
only marginally distinctive. 

The second set of panels shows the azimuthal distance between the leading tagging jet and the heavy res- 
onance. The Standard Model Higgs and the cut-off 2g W structures predict the same moderate back-to-back 



topology. For spin-1 the l v 
kinematics. 



and 1Z, distributions are much more strongly peaked in the same back-to-back 



Comparing the jet- A correlations shown in Fig. [9] with the jet-jet correlations from Fig. [8] clearly shows that 
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Figure 9: Normalized correlations between the leading tagging jet and the heavy resonance. We show the difference in 
the rapidity difference Arjjx as well as the azimuthal angle A(j>jx- The spin-0, spin-1, spin-2 interaction structures are 
described in Sec. IIIIAIIIICI 
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in weak boson fusion the spin and coupling information is largely encoded in the tagging jets, not in the angular 
correlations of the heavy resonance [TT], with the notable exception of a unitarized spin-2 model [T2"] . 



C. Including Higgs decays 



Essentially all information on the Higgs couplings structure available from 'Higgs' production in weak boson 
fusion is included in the jet-jet and jet-X correlations described in the previous two sections. The only caveat 
is that the resonance observed in the final state might have a spin or polarization structure which we average 
over unless we make certain requirements on the X — > dd decay products. The price we have to pay for this 
additional information is a full simulation of the X decay, including its underlying coupling structure. If we 
limit ourselves to VVX couplings for the production process this obviously implies additional assumptions for 
all decays except for decays to massive weak gauge bosons X — > VV. Moreover, the production and decay 
process VV — > X — > VV requires a unitarization which only the 0g M delivers and which we have to add 
for example for a spin-2 resonance. In Fig. [7] we see that this unitarization has a major effect on all energy 
dependent distributions, which means we have to take the jet-decay correlations with a grain of salt. 



Nevertheless, in Fig 10 we show the usual rapidity and azimuthal angle differences between the leading 
tagging jet and one of the two X decay products. For the decay channel we assume X — > t + t~ with the 
couplings described in Sec. |III| Comparing these distributions to the jet-JT distributions in Fig. [9] we see that 
their features are very similar, i.e. the heavy resonance decay kinematics is dominated by the jet-X kinematics. 
Small effects come from the acceptance cuts for the decay products. The main target for production-decay 



correlations is the spin-2 hypothesis, which after the unitarization cut pxj < 
dominated by the helicity-2 states. 



..max 
FT 



in Eq.(20) tends to be 




Figure 10: Normalized correlations between the leading tagging jet and the X decay particle. We show the difference 
in the rapidity difference Ar/jd as well as the azimuthal angle A(/>jd- In the bottom row we show the angle <£>+ in the 
Breit frame. 
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The spin-2 operator coupling to gauge bosons and to quarks is the only scenario where Arjjd differs from Ar/jx- 
While the latter shows a centrally flat behavior, distinctively different from 2^ w and 0g M , the distributions 
including the decay becomes much more similar to those. In this case the jet-decay correlations are even less 
useful than the jet- A correlations. The distributions also resemble their A<f>jx counter parts closely, 

with a generic smearing because the tagging jets now recoil against to decay products instead of one resonance. 
Again, one change is that the two spin-2 hypotheses are becoming indistinguishable, but possibly more different 
from the spin-0 predictions. Moreover, one exception to the general broadening of A(f>jd is the \\ operator, 
which becomes more similar to the other spin-1 scenarios and less similar to the scalar couplings. 

In addition to the generic hadron collider observables there are a few variables which are particularly well 
suited to identify a spin-2 resonance. One example is the Gottfried-Jackson angle [31] [32], which is very closely 
related to cos 9* in the Breit frame. We show the corresponding distributions in the Appendix. In the bottom 



row of Fig. 10 we show the distribution of the angle $+ = 2<j>i + A0 in the Breit frame [12] . whose angles are 
defined in Eq.Q. Its particular feature is that it has no distinguishing features unless we unitarize the spin-2 
rate, in which case it develops a clear modulation for spin-2 couplings. In principle the same information should 
be included in an appropriate combination of hadron collider observables, but it cannot be easily observed in 
combinations like A(f>j 1( i + A<j)j 2 d- 



D. Basic strategy 

From the discussion of the different jet-jet, jet-A, and jet-decay observables we know that there are a few key 
measurements which can distinguish between the 'Higgs' coupling structures. The hadron collider observables 
Ar] and A(j> serve as effective replacements of the Cabibbo-Maksymowicz-Dell'Aquila-Nelson angles [131 114) . 
A particular strength of the angular basis for X — > ZZ decays is the clear link between some of some Higgs 
couplings structures and particular distinctive observables. The distribution of the azimuthal angle A<f>jj |llj 
shown in Fig. [2] constitutes such a specific link in the WBF topology. 

Hence, the question arises if we can construct a set of measurements which can distinguish the different 
'Higgs' operators shown in Tab.[lj Such a simple flow diagram does not correspond to an optimized experimental 
analysis, but it indicates what the key observables are. For example, from Sec. |IV| we know that we should avoid 
cutting on Arjjj before testing 'Higgs' couplings because it removes one of the key observables to tell apart the 
0p 5 and 0g M operators. Following our earlier argument we do not include the transverse momentum spectrum 
of the tagging jets in our list of observables because the initially striking differences will be modified by form 
factors for spin-1 and spin-2 resonances. For now, we consistently assume spin-2 operators with our simplified 



recognized from the pt spectra shown in Sec IV A 



militarization cut p™ ax = 100 GeV as defined in Eq.(20). Without this cut a spin-2 could be immediately 



In Fig. [TTJwe show some of the most distinctive measurements reflecting the couplings of a heavy resonance. 
The first jet-jet variable, A<f)jj, is well suited to clearly separate the spin-0 structures. As we will see in the 
following section the Arjjj distribution might be even more powerful than Atfijj , so there should be no problem 
in distinguishing the different spin-0 structures just based on the tagging jets. 

The problem is to tell apart some of the spin-1 and spin-2 hypotheses from their closest spin-0 models in the 
jet-jet correlations. For example, to tell apart the Standard Model Higgs boson 0g M and the graviton 2g W+q 
we have to resort to the jet-A correlations Arjjx- From Sec. IV C we know that Breit frame angles would also 
be well suited for this distinction, most notably <fr+ [12) . 

Similarly, telling apart 0^ 5 and \\ we requires jet-decay correlations like Arjjd- Alternatively, we could use 



4>+ shown in Fig. 10 the Breit-frame angle 9*, or the Gottfried- Jackson angle 9qj. The two operators l w z 
appear indistinguishable using the set of observables at hand. 

The qualitative bottom line of this section is that clearly the jet-jet correlations are the most decisive, in 
particular to separate the different scalar coupling structures [11) . For some spin-1 and (unitarized) spin-2 
models we need to add information from the heavy resonance or its decays [T2] . 



E. Comparison of observables 



Following Sec. IV D we know what the candidate observables for the distinction of model hypotheses are. 



However, for example for the distinction of the + 



SM 



and 0p 5 coupling structures there are several such variables. 



Only a complete shape analysis can tell us how promising these distributions really are. 
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The statistical power of a well-defined set of observables we compute based on statistical errors, ignoring 
systematics as well a theoretical uncertainties and errors. This approximation is justified for relatively low 
luminosities — for example asking the question which of the distributions would first cross a given confidence 
level in the comparison between two hypotheses. 

To asses such a minimal luminosity we perform a binned log-likelihood ratio hypothesis test |47j . Our zero 
hypothesis or background hypothesis is the Standard Model coupling structure 0g M . This hypothesis also fixes 
the number of events in the distribution. The exotic signal we would like to discover or reject are the alternative 
spin-0 or spin-2 coupling structures. Backgrounds we ignore in this computation, because they mostly affect 
the exact value of the luminosity required for a 95% exclusion, but not the relative strength of the different 
observables. The binned likelihood is defined in terms of a Poisson likelihood L, 

N n (X)e- N( - x 1 
Lfdatal hypothesis X) = ; 

TV. 

n- ok. L ( data I hypothesis 1) 
y_ g L(data|SM) ' 

The expected number of signal events (per bin) given the luminosity C is N(X) — a(X)C, while n is the observed 



number of signal events. The logarithm is additive, so the generalization of Eq.(21) to binned distributions 



with (in our case) 20 bins per distribution is straightforward. The Neyman-Pearson lemma states that this is 



the statistically 'best' discriminator. The different hypotheses in Eq.(21 ) are sampled with randomly generated 
pseudo-data, giving us probability density functions whose overlap define the 95% CL curves. As alluded to 
above, the absolute value of the luminosity required for a 95% CL exclusion of alternative coupling structure 
should be taken with a grain of salt. What we aim at is the relative strength of the different observables, which 
might be hard to judge by eye. 

In Fig. [12] we compare the statistical power of different observables in WBF kinematics. In the left panel we 
show how well the different observables can distinguish the 0p 5 operator from the Standard Model coupling 
structure. The two most promising observables are the rapidity difference between the two tagging jets and 
their azimuthal angle difference, i.e. observables which are independent of the Higgs decay. All that is required 
is that the event sample has a good signal-to-background ratio, such that such an analysis is not dominated 
by the systematics of the background subtraction. The necessary clean-up of the event sample can be achieved 
by central jet vetos [37] and/or the matrix element method [3U]. We see that the rapidity difference turns out 
more powerful than the azimuthal angle. 

In the second panel we show the same set of observables for a discovery of the CP-odd 0p 5 coupling structure, a 
pseudo-scalar Higgs boson. Again, the azimuthal angle A<pjj is somewhat less sensitive compared to the rapidity 
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Figure 12: Confidence level for the distinction of different spin-0 coupling structures as a function of the integrated 
luminosity. We assume statistical uncertainties only for the WBF signal in the absence of backgrounds. The plot- 
ted confidence level refers to exclusion of the SM WBF component. The cross section normalization is 2.8 fb which 
corresponds to a good r reconstruction in the purely leptonic channels (see |19|). 



difference. Finally, we study the most relevant observables to identify a 2 + resonance without a constraint on 
the high-energy tails of the jet momenta. As before, the rapidity difference of the tagging jets is the most 
powerful. These results are a direct consequence of the high discriminative power of the Arjjj distribution 
shown in Fig. [8j where all the alternative couplings tested in this section peak around Arjjj ~ 4, compared to 
the Standard Model 0g M case with Arjjj ~ 5.5. For spin-2 identification adding information from the heavy 
resonance or its decays implies a significant improvement |12) . Our results for the best-suited distribution 
reflects the setup of the likelihood test: the required luminosity is dominated by those bins which show a large 
ratio of events between the two different hypotheses. 

Going beyond our idealized setup, these jet-jet observables differ when it comes to detector and QCD effects, 
e.g. pile-up suppression is a function of rapidity Furthermore it is well-known that A<f>jj acts as good 
measurement quantity to extract mixed CP properties [2U] • Such an angular correlation is insufficiently reflected 
in the rapidity difference. This way Aipjj remains a well-suited observable to test more complex CP properties 
in WBF Higgs production, while a wide operator spectroscopy preferably relies on the leading observable Aijjj . 



V. ZH PRODUCTION 



The study of WBF kinematics with the different decay channels gives us a wide array of tests of the Higgs 
coupling structures. However, fermions still only feature in one channel, namely X —> t + t~ decays. Once the 
LHC runs at 14 TeV, another channel will become available to test the 'Higgs' couplings to fermions, 



PP 



(Z^£+£-)(X -> bb) 



(22) 



Independent of the use of modern Higgs finding methods [21] this channel will only be observable for boosted 
Z and X decays. The kinematics of this process can be directly mapped onto the decay X — > ZZ, so the 
Cabibbo-Maksymowicz-Dell'Aquila-Nelson angles [13] are a natural choice for spin-sensitive observables. In 
complete analogy to Eq.Q we define 



PZ' =Pz +Px 



Px =Pb+ Pi 



Pz 



e- +P?+ 



(23) 



as well as the corresponding unit three-momenta pi in the Z* and the Z, X rest frames. Similar to the original 
angles defined in Sec. |II A] the Z* in the s-channel will be off-shell, which does not prevent us from using it as 
a reference frame. The basis of angles is then defined as 



cos u b = p bl ■ p z 



x 



PX ■ Pbcam 



cos 



4>b = (Pbcam X pz) ■ (pz X p bl ) 



COS Vi = Pi- ■ p X 

cos = (p b x p- b ) ■ (pg- x p l+ ) 

The signal distributions for the different X hypotheses are generated with the rough acceptance cuts 
PT,b > 20 GeV p T ,£ > 10 GeV AR > 0.4 |r? 6 | < 5 |»#| < 2.5 



(24) 



(25) 




where we assume the same value for all AR value for 6-jets and leptons. In addition, we require 

p Tm > 100 GeV , (26) 
to be consistent with the proposed LHC analyses. 

Because this analysis requires reconstructed Z — > £ + £~ and the H — > bb decays we limit our model hypotheses 
from Tab.|l]to those which predict a XZZ coupling as well as a Xbb coupling. For this sub-set of spin-0, spin-1, 



and spin-2 models we should the most discriminating distributions in Fig. 13 The three scalar operators can 
be identified by a flat cos 9b distribution. The SM and the two higher-dimensional scalar couplings then show a 
distinctively different behavior in cos 9i, which carries the same physical information as the cos# e distribution 
presented in Fig(3] for the X — > ZZ decay. The challenge in the scalar sector will be the measurement of the 
CP property of the 0^ 5 operators. The problem with these two states is that they are both scalars, washing 
out any angular correlation between the X — > bb and the Z — > £ + £~ decays. If their dimension-5 origin gives 
them a similar energy behavior we would need a dedicated CP asymmetry to distinguish them. 

Comparing the different spin-1 and spin-2 hypotheses is made easy by the distinctive forward-backward peaks 
in the spin-2 prediction for cos^. For 9i and 8* the 1^ distribution has a significantly larger preference for 
cos 9 = 0. While the boosted kinematics might well start to wash out the angular correlations of the Z, X 
decay products, this simple set of predictions indicates that the determination of Higgs coupling structures in 
associated ZH production is promising. 



VI. OUTLOOK 



After the discovery of a Higgs-like resonance by ATLAS and CMS the main focus of Higgs analyses at the 
LHC will be the determination of the Higgs Lagrangian. This includes the structure of the operators (linked 
to the spin and CP quantum numbers of the 'Higgs' boson) as well as an independent measurement of the 
coupling strength. Similar to the hugely successful electroweak precision program centered around LEP this 
should give us hints about the embedding of the Higgs mechanism in an ultraviolet completion of the Standard 
Model. 

We present a comprehensive study for the determination of the Higgs coupling structure in weak boson 
fusion. Starting from the Cabibbo-Maksymowicz-Dell'Aquila-Nelson angles in X — > ZZ decays we show how 
the WBF kinematics can be described either using Breit frame angles or rapidity and azimuthal angles. 

The most distinctive observables in the WBF topology are jet-jet correlations. It is well known that the 
azimuthal angle A<fijj can distinguish between different spin-0 couplings without ever requiring a reconstructed 
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resonance momentum |1U 120) . Based on a signal-only and parton level statistical analysis we expect the 
rapidity difference Arjjj to be at least as effective in distinguishing different spin-0 and spin-1 hypotheses. The 
advantage of jet-jet correlations is that we do not have to model 'Higgs' decays based on different models, where 
the production in weak boson fusion and a decay to photons or fermions might probe very different underlying 
assumptions on the production and decay sides. 

Generalizing our analysis to spin-1 and spin-2 hypotheses requires additional information from the heavy 
resonance decays |12j . Again, the more general approach is to reconstruct the heavy resonance momentum 
without any assumption on the decay vertex and test angular jet- A correlations. We find that they indeed 
allow for a distinction of spin-0 and spin-1 or spin-2 models. Adding decays of the heavy resonances and 
computing jet-decay correlations adds only little in terms of the hadron collider observables Arj and Acf>. 

Nevertheless, it might be useful to add dedicated angular correlations defined differently. This includes an 
analytically continuation of the X — > ZZ angles, or some Breit frame angles like 9* or <3? + . 

For associated ZH production with fully reconstructed decays Z — > £ + £~ and X — > bb we can immediately 
employ the X — > ZZ angles. It should be possible to confirm the 0g M coupling structure also in this channel 
without any problems. 
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Complete set of observables 



In Sections IV A IV B[ and |IV C" we show a selected set of observables which should allow us to determine the 



coupling structure of a heavy Higgs-like resonance in weak boson fusion. We largely rely on the hadron collider 
observables Aij and A(j>, even though in Sec jll A| we introduce Breit frame observables as well as dedicated 
angles like the Gottfried- Jackson angle. In this appendix we give a complete set of distributions for all model 
hypotheses listed in Tab[TJ 

In Fig. [14] we start by giving the complete set of rapidity and transverse momentum distributions for the 
heavy resonance X and the leading tagging jet. In addition to the curves shown in Fig. [7] we also show the 
scalar and spin-2 couplings to incoming gluons, as defined in Tab. [I] 



spin-0 spin-1 spin-2 spin-2(p™ ax ) 



O.S r I 0.5 r ] 0.5 r i 0.5^^ 




Figure 14: Normalized distributions of the heavy resonance X and the leading tagging jet after the minimal kinematical 
cuts of Eq.(17l and rrijj > 600 GeV in Eq.(18l. The hypotheses are 0g M (red), 0q 5 (blue), 0q 5 (green), 0^ 5 , , (blue 



dashed), D5 j g ^ (green dashed); lj (blue), (blue dashed), l z (green), l w (green dashed); 2^ w (black), 2 



dashed), 2q CD (blue), 2 + (blue dashed), all defined in Tab. |IJ 



D5(g) 
EW+q 



(black 
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spin-0 spin-1 spin-2 spin-2(p™ ax ) 




Figure 15: Correlations between the tagging jets. The different curves are the same as in Fig. [14] 



In Figs. 15J THJ and 17 we show the jet-jet, jet-X, and jet-decay correlations for the same full set of operators. 
The confirm the picture that comes out of Sec. |IV[ that for the WBF topology the jet-jet correlations are the 
most sensitive because they show the most distinctive kinematic patterns. Moreover, Fig. [15] confirms that 
adding a small fraction of gluon fusion events hardly changes or even dilutes the general features of the WBF 
kinematics. For all distributions we also observe a qualitative difference between the generic spin-2 couplings 
and the spin-2 couplings after a unitarization cut px,j — p™ ax . A completely model independent test of spin-2 
couplings is unfortunately not possible. 



Finally, in Fig. [18] we show the most promising specific angular observable for the determination of the 
coupling structure in weak boson fusion. This includes the Breit-frame angles $+ and cos 9* defined in Eq.Q 
as well as the Gottfried-Jackson angles introduced at the end of Sec. |II A] We see that their benefit apart from 
a spin-2 identification is somewhat limited. 
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Figure 16: Same as Fig. |15[ but among the heavy resonance X and the leading tagging jet. 




Figure 18: Specific angles targeted at identifying a spin-2 resonance. The different curves are the same as in Fig. |14| 
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